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ABSTRACT 

Very low-mass binaries (VLMBs), with system masses < 0.2 M appear to have very 
different properties to stellar binaries. This has led to the suggestion that VLMBs form 
a distinct and different population. As most stars are born in clusters, dynamical evo- 
lution can significantly alter any initial binary population, preferentially destroying 
wide binaries. In this paper we examine the dynamical evolution of initially differ- 
ent VLMB distributions in clusters to investigate how different the initial and final 
distributions can be. 

We find that the majority of the observed VLMB systems, which have separations 
< 20 au, cannot be destroyed in even the densest clusters. Therefore, the distribution 
of VLMBs with separations < 20 au now must have been the birth population (al- 
though we note that the observations of this population may be very incomplete). 
Most VLMBs with separations > 100 au can be destroyed in high-density clusters, 
but are mainly unaffected in low-density clusters. Therefore, the initial VLMB pop- 
ulation must contain many more binaries with these separations than now, or such 
systems must be made by capture during cluster dissolution. M-dwarf binaries are pro- 
cessed in the same way as VLMBs and so the difference in the current field populations 
cither points to fundamentally different birth populations, or significant observational 
incompleteness in one or both samples. 
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1 INTRODUCTION 

It has been suggested that brown dwarf-brown dwarf bina- 
ries, or, more generally, very low mass binaries (VLMBs) 
with system masses < 0.2 Mq form in a different way to 
stellar binaries. The main argument for this scenario is that 
the binary fraction and separation distributions of VLMBs 
are very different to those of stars. iThies fc" Kroupal (|2007n 
point out that the binary fraction of very-low mass sys- 
tems is only 15 - 25 per cent, compared to 42 per cent for 
higher-mass M-dwarfs. Also, the separation distribution for 
most stellar binaries (higher-mass M-, K-, and G-dwarfs) 
has the same mean (30 au) and variance (o"i oglo a = 1.53, 
where a is the semi-major axis in au), and differs only 
in the multiplicity of the primary in the particular mas s 
range ijDuauennov fc Mavorl [l99ll : iFischer fc Marcvlll992h . 
However, the observed separat ion distribution of VLMBs 
(see e.g. iBurgasser et al.l 120071 ) show the data (when fit- 
ted with a logio-normal) to have a mean of 4.6 au with a 



E-mail: rparker@phys.ethz.ch 



much smaller variance. In addition, IThies fc Kroupal (|2007l . 
2008) argue that the observations of VLMBs are not con- 
sistent wotli_a_contnmous^ over the hydrogen-burning 
limit. IThies fc Kroupal J2007) interpreted this as evidence 
that VLMBs form through a different mechanism to stellar 
binaries. 

However, it is known that binary populations can un- 
dergo significant dynamical processi ng, with many, espe- 

"" 19751: iKroupal 
Kroupa et al.l 



cially wi der systems, being destroyed dHeggid 



| l995a|lbl:lKroupa. Petr fc McCaughreanll 19991 : 



20031 : [Parker et al.ll2009l ). Therefore, the currently observed 



binary population, es pecially in the fi eld, is not the same as 
the birth 



npulation, especiall y m the f ie 
population (|GoodwirJl20"l~bT ). 



In this paper, we examine to what extent an initial 
VLMB population can be altered by dynamical process- 
ing and so how different the initial VLMB and stellar bi- 
nary populations ca n be at birth. We extend the work of 
iKroupa et all (|2003l ) who examine the evolution of a mixed 
population of star and very low-mass object (VLMO) bina- 
ries and find that, dynamicall y at least, VLMBs must form a 
separate population (see also lThies fc Kroupall2007l [2008). 
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Here we assume that VLMBs are a separate population and 
examine what range of initial binary fractions and separa- 
tions can reproduce the current observations. In Section [2] 
we review the available VLMB data, in Section[3]we describe 
the set-up of our simulations, we present our results and dis- 
cussion in Sections [4] and [5] and we conclude in Section [6] 



2 SUMMARY OF THE OBSERVATIONS 

Data on very low mass binaries (VLMBs) - binaries with 
a total system mass of < 0.2 Mq - have been collated 
in the Very Low M ass Binaries Archive (VLMBA, see 
iBureasser etafl [20070 Given that the majority of these 
systems have primary masses less than the hydrogen- 
burning limit (m H = 0.08 M ), the VLMBA data have the 
potential to provide an excellent constraint on the hypoth- 
esis that substellar binaries form via a different mechanism 
or in a different environment to stellar binaries. 

As of September 2010, the VLMBA lists 99 systems, 
four of which lack robust separation measurements and a 
further system has a planetary-mass companion. That leaves 
94 systems to compare with numerical simulations. 



2.1 Multiplicity 

We define the multiplicity of VLMBs as 
B 



/VLMB = S + B , (!) 

where B is the number of binary systems and S is the num- 
ber of single systems. We ignore triple and higher-order sys- 
tems for the remainder of this paper. 

The overall multiplicity of VLMBs is open to debate. 
Based on potentially undiscovered systems. iBasri fc Reinersl 
(2006) suggest a a value of 0.26 ± 0.10. This would argue in 
favour of the VLMBs as a continuous population, as the mul- 
tiplicit y of M-dwarfs and G-dwarfs is 0.42 and 0.58, respec - 
tively (|Fischer fc Marcvl 1 19921 ; iDuquennov fc Mavorl Il99lh 
possibly indicating a smooth decre ase in multiplicity with 
decreasing primary mass. However, iThies fc Kroupal \200H ) 
argue that the VLMBA data is consistent with an overall 
multiplicity of 0.15, representing a distinct cut-off from the 
stellar binary regime. 



2.2 Separation distribution 

In Fig. [T] we plot the separation distribution of the observed 
VLMBs, distinguishing between systems observed in the 
Galactic field (the hashed histogram) and systems observed 
in various clusters (the open histogram). The majority of 
these systems lie in a small separation range (~ 1 — 30 au 
with a peak around 5au). The total area of the histogram 
corre sponds to the obse rved binary fraction of VLM objects, 
0.15 l|Close et alj|2003l ). 

It is important to note that most of the observed 
VLMBs are in the field. This means that their birth cluster 
has disrupted and that we are mostly observing an already 



1 See http:/ /www.vlmbinaries.org/| for an up-to-date census of 
the known very low mass binaries, maintained by N. Siegler, 
C. Gelino and A. Burgasser. 



Table 1. Parameters for the logio-normal fits to the data in 
Fig. [TJ The columns contain, from left to right, the reference, 
mean, log of the mean, variance and the binary fraction used to 
norma lise the distribution. The refere nces arejThie^^^Krouga 



J2007I . TK07), IBasri fc Reinersl d2006l. BR06) , iFischer fc Marcv 
l|l99a . FM92) and lDuauennov fc Mavorl il99ll . DM91). 



Ref. 


Mean 


Mean 


Variance 


/mult 




(au) 


(logio a) 


(°log 10 a) 




TK07 


4.6 


0.66 


0.40 


0.15 


BR06 


4.6 


0.66 


0.85 


0.26 


FM92 


30 


1.57 


1.53 


0.42 


DM91 


30 


1.57 


1.53 


0.58 



dynamically processed VLMB population. Unfortunately we 
have no information on the birth clusters of these systems. 

Several authors have attempted to fit the VLMB 
data with logio-normal distributions similar to those 
that a re used to fit G-, K- a n d M-dwarfs in the 
(IDuquennov fc Mavorl Il99ll ; iMavor et ail 

fir 



field 



I1992I; 



IFischer fc Marcvl 1 199^ respectively ]j. IThies fc Kroupal 
it the VLMB separation distribution with a 
logio-normal distribution with mean 4.6 au 



O"log 10 a 



and variance 
0.4, and an overall substellar multiplicity of 0.15 



(the (solid) brown line in Fig. [TJ . 

However, there are some outlying systems with very 
short and long separations, and lBasri fc Reinersl (|2006t ) ar- 
gue for a wider distribution, based on the hypothesis that 
there may be unresol ved short-period VLMBs with separa- 
tions less than 1 au (|Maxted fc Jeffries! I2005T ) and (at the 
time) the tentative d iscovery of VLMBs with separations in 
excess of 100 au (e.g. IClose etaill2003l ; iBouv et alj|2006j ). 

The former is still the subject of debate, with claims 
that the peak in t he VLMB separation d istributio n may be 
betwe en 1 - 3 au l|Burgasser et alj l2007h , though Ijoergensl 
(2008) suggests that few VLMBs exist with separations 
<3au. The latter appears to be partly vindicated by re- 
cent discoveries o f wider systems in the field (Konigstuhl-l 
AB, a = 1800 au.lCaballeroll2007l ; 2M0126A B a = 5100 au 
Artigau et alj|2007l : 2M1258AB. a = 6700 au. lRadigan et all 
2009), although surveys should b e sensitive t o VLMBs with 



sep arations between 1 - 2 00 au (|Burgasser et al.ll2007h . 

IBasri fc Reinersl (|2006l ) proposed a wider logio-normal 
fit to the data with mean 4.6 au and variance o"i O g 10 a = 
0.85, and an overall substellar multiplicity of 0.26 (the 
(dot-dashed) magenta line in Fig. [TJ. For comparison, in 
Fig. [TJ we also show the logio-normal fits for field M-dwarfs 
(IFischer fc Marcvll 19921. the (dashed ) blue line) and field G- 
dwarfs (|Duquennov fc Mavorl Il99ll . the (dotted) red line). 
Details of the parameters for the logio-normal fits are given 
in Table [TJ 



2 Whilst the fit to the G-dwarf data is logio-normal, the K- and 
M-dwarf separation distributions suffer from poorer statistics and 
may not be a straightforward scaling down of the G-dwarf distri- 
bution as a function of the multiplicity of the primary component 
(as multiplicity decreases as a function of primary mass in bina- 
ries in the field). 
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Figure 1. Data from the Very Low Mass Binary Archive 
(VLMBA - iBurgasser et alj 120071 ). The hashed histo gram rep- 
resents binaries observed in the Galactic field; whereas the open 
histogram represents the few VLMBs observed in various star 
clusters. The total area of the histogram corresponds to the ob- 
served binary f raction of VLMBs, 0.1 5. The logio-normal fits 
to t he data by iThies fc Kroupal ||2007| , the (solid) brown line) 
and iBasri fc Reinerj ( 1200^7 the (dot-dashed) magenta line) are 
shown. For comparison , the logio-normal fits for field M-dwarfs 



Fischer 



Marcvll992l. the (dashed) blue line) and field G-dwarfs 
Duquennov fe M ayor 1991, the (dotted) red line) are also shown. 




Figure 2. The mass ratio distributio n of systems from the Very 
Low Mass Binary Archive (VLMBA - IBurgasser et alj|2007l) . The 
bins are of width 0.1, and are normalised to the total number of 
systems used (94) to allow direct comparison with the simula- 
tions. 



2.3 Mass ratio distribution 

It is also interesting to trace the possible evolution of the 
mass ratio distribution. For each system, the mass ratio, q, 
is denned as 

g=^, (2) 

where m p and m s are the masses of the primary and sec- 
ondary components, respectively. 

In Fig. we show the observed mass ratio distribution 
of the VLMBA data, normalised to the total number of sys- 
tems (94). Almost half the VLMBs in the sample have a 
mass ratio approaching unity, and the majority of the other 
systems have high (> 0.7) values of q. 



2.4 Other properties 

Data on other dynamical properties of the VLMBs are not 
yet available. Therefore, we do not include a study on e.g. 
the possible effects of cluster evolution on the VLMB eccen- 
tricity distribution. 



3 METHOD 

3.1 Cluster set-up 

We follow a simi lar method to the one described in 
IParker et all l|2009D to set up the clusters and stellar bi- 
nanjf] systems in our simulations. The clusters are designed 
to mimic a 'typical' star cluster, similar to Orion with 
N = 2000 members and mass ~ 10 3 M Q . 

For each set of initial conditions, we create a suite of 
10 simulations, corresponding to 10 clusters, identical apart 
from the random number seed used to initialise the simula- 
tions. 

We set our clusters up as initially virialised 
Plummer spheres JplummeJ 191 if ) as described by 
l Aarseth. Henon fc Wielenl | | 19 74 )7" The prescription in 
lAarseth et al.l 1 19741 ) provides the positions and velocities of 
the centres of mass of the systems in the Plummer sphere. 
The current half-mass radius of Orion is 0.8 pc 



1 McCaughrcan & Stauffei] fl994l; lllillenbrand fc Hartmannl 




Kohler et alj 120061 ) 



argue that Orion was 



However, 
originally 



Parker et al.l 



much denser 



sion (iTutukovl 1 19781; 


Hills! Il98d; iGoodwinl Il997l; 


Krouoa, Aarseth & Hurley 


2001; iGoodwin & Bastian 


2006) and dynamical 


interactions JKrouoa et al. 


19991, Ivan den Berk, Porteries Zwart & McMillanl 20071; 


Parker et al.l 20091, and references therein) have caused it 



to expand to its current size. We therefore adopt initial 
half-mass radii of 0.1 pc and 0.8 pc for the clusters in 
our simulations, thereby covering a wide range of cluster 
densities. 



3 From this point in the paper, we adopt the phrase "stellar bi- 
nary" when describing systems with component masses both ex- 
ceeding 0.106 Mq; and "very low-mass binary (VLMB)" when de- 
scribing systems with both component masses less than 0.106 Mq. 
0.105 Mq is the mass of the most massive VLMB primary com- 
ponent in the archive. 
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Observations suggest that the ra tio of stars with masse s 
< 1 Mq to brown dwarfs is ~5:1 (e.g. lAndersen et al.| [2008). 
In our simulations, we place one sub-stellar system (either 
single or binary) in the cluster for every five stellar systems. 

3.2 Stellar binary properties 

It is thought that the star formation pr ocess should produce 
binary stars in preference to singles l| Goodwin fe Kroupal 
120051 ; [Goodwin et al.ll2007l . and references therein). There- 
fore, all the clusters in our simulations are formed with an 
initial stellar binary fraction, / ate iiar = 1 (i.e. all stars form 
in binary systems; there are no singles or triples, etc.), where 

B 



e tid = - [0.95 + tanh (0.6 logioP - 1.7)] . 



(7) 



S + B' 



(3) 



and S and B are the numbers of single and binary systems, 
respectively. 

The mass o f the primary star is chosen randomly from 
a lKroupal (|2002T ) IMF of the form 



N(M) oc 



M~ 13 mi < M/M < m 2 , 
M~ 2 3 m 2 < M/M < m 3 , 



(4) 



where mi = 0.106 M , m 2 = 0.5 M , and m 3 = 50 Mq. 
Note that the lower mass limit, mi, is higher than in our 
previ ous papers (e.g. iParker fc Goodwinl l2009; Parker et al.l 
120091 ). This is to prevent a stellar binary or single star from 
having mass components (m p and m s ) that would overlap 
with the VLMBA datfl 

In several test simulations we found 
that varying mi by a few per cent causes a negligible differ- 
ence to the results. 

Secondary masses are drawn from a flat mass ratio dis- 
tribution with the constraint that if the companion mass is 
< 0.106 M or the total mass of the binary, m to t < 0.2 M 
it is reselected, thereby removing the possibility of creating 
VLMBs in the stellar domain. Note that this method does 
not produce the input IMF, exactly due to the way that 
secondaries are chosen. 

In accordance wi th the observation s of 

iDucmennov fc Mavorl (|l99lh and iFischer fc Marcvl <|l992h . 
the periods of stellar binary systems are drawn from a 
logio-normal distribution of the form 



/ (logioP) oc exp 



- (logio P - logio P) 

2 °"log 10 P 



(5) 



where logioP = 4.8, <7i oglo p = 2.3 and P is in days. The 
periods are then converted to semi-major axes. 

Eccentricities of stellar bi naries are drawn from a ther- 
mal eccentricity distribution (|Heggidfl975l ; IKroupal [l995al , 
2008) of the form 



Me) = 2e. 



(6) 



Binaries with small periods but large eccentricities 
would expect to undergo the tidal circularisation show n in 
the sample of G-dwarfs in iDuquennov fc Mavorl (| 199 if ) . We 
account for this by reselecting the eccentricity if it exceeds 
the following period-dependent value etid: 



4 We note that in their work, IThies fc Kroupal | |2007|) deliber- 
ately allowed an overlap in the mass range 0.08 — 0.1 M , so that 
objects could be either stars or brown dwarfs. 



with logio P in dayfl This ensures that the eccentricity- 
period distribution ma tches the observations of 
IDucmennov fc Mavorl (|l99ll ). 

3.3 Very low-mass binary properties 

In three sets of sim ulations, we use the VLMBA data 
( Burgasscr et al. 200^) to randomly choose binaries and use 
their masses and semi-major axes as initial values for sub- 
stellar binaries in the cluster. 

We run simulations wit h semi-major axes drawn from 
the IThies fc Kroupal (120071 ') fit to the VLMB separation 
distribution, with Gaussian parameters of logio a = 0.66, 
where a is the semi-major axis in au (0.66 corresponds 
to 4.6 au), and var i ance o"i ogln a — 0.4; and also the fit 
by iBasri fc Reinersl (|2006T ). which accounts for the outly- 
ing binaries in the VLMB separation distribution by adopt- 
ing the same logio-normal peak but increasing the variance 
(logio a = 0.66, a\ ogloa = 0.85). For completeness, we run 
simulations in which the sub-stellar binaries have the same 
separation distribution as the stellar binaries (logio P = 4.8, 
fio gl0 p = 2.3; P in days). 

In each case, we adopt an initial VLMB fraction; either 
0.5 (c .f. the stellar binaries in the field lDuauennov fc Mavorl 
Il991 ), 0.25 (to accoun t for potential ly undiscovered VLMBs , 
IBasri fc Reinersll2006r ). or 0.15 (the lThies fc Kroupal (|2007h 
fit to the observations). 

In the simulations that choose separations from the vari- 
ous logio-normal distributions, the masses of the sub-stellar 
binary components are chosen by randomly assigning the 
primary a mass in the range 0.01 M < m p ^ 0.106 M . We 
then adopt a flat mass ratio (q) distribution to choose the 
mass of the secondary component (0.01 M ^ m s < m v ). 

This mass range allows a direct comparison between 
the VLMBA and the simulations. We find that reducing 
the upper limit of the VLMB primaries to 0.08 M has a 
negligible effect on the results. 

The eccentricities are drawn from a thermal eccentricity 
distribution and then tidally circularised (Eqns.[6]and[7]). 

A summary of the different mass, semi-major axes and 
eccentricity distributions used to create the substellar bina- 
ries is given in Table [2] 



3.4 V-body integration 

By combining the primary and secondary masses of the 
stellar and VLMBs with their semi-major axes and eccen- 
tricities, the relative velocity and radial components of the 
stars/very low-mass objects in each system are determined. 
These are then placed at the centre of mass and centre of 
velocity for each system in the Plummer sphere. 

Simul ations are run using the kira int egrator in Star- 
lab (e.g. iPortegies Zwart et al l Il999l . l200ll . and references 



5 A more elaborate 'eigenevolution' algorithm, which accounts 
for pre-Main Sequence ti dal circ ularisation and protostcllar disk 
accretion is described in IKroupal j2008h . We elect not to use it 
here, as the disk accretion mechanism in this algorithm also alters 
the mass ratio distribution. 
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Table 2. A summary of the different substellar binary con- 
figurations adopted in the simulations. Prom left to right, the 
columns show the initial cluster half-mass radius, ri/2'1 the 
se paration distribution (either the log ip-normal distributions 
of iDuquennov fc Mayor] il99ll. DM 91), iBasri fc Reinersl ||2006| . 



BR06) and iThies fc Kroupal J2OO7I , TK07) or values taken di- 



rectly from the VLMBA); the initial multiplicity; the masses of 
the two components in a system (either randomly selected with 
a flat mass-ratio distribution, or taken from the VLMBA). See 
Section 13.31 for full details. 



ri/2 


Separation Distribution 


/VLMB 


Mass Range 


0.1 pc 


VLMBA 


0.25 


VLMBA 


0.8 pc 


VLMBA 


0.25 


VLMBA 


0.8 pc 


VLMBA 


0.15 


VLMBA 


0.1 pc 


TK07 


0.15 


Random, flat q 


0.1 pc 


BR06 


0.5 


Random, flat q 


0.1 pc 


BR06 


0.25 


Random, flat q 


0.1 pc 


DM91 


0.5 


Random, flat q 


0.8 pc 


DM91 


0.5 


Random, flat q 



therein) and evolved for lOMyr. We do not include stellar 
evolution in the simulations. 



4 RESULTS 

4.1 Finding bound binary systems 

We use the neares t-neighbour algorithm described by 
Parker et alj d2009h (a nd independently verified by 
Ll[2( 



Kouwenhoven et al.ll2010l ) to determine whether a star/very 
low-mass object is in a bound binary system. 



4.2 The evolution of the VLMB separation 
distributions 

In this section we discuss the effects of cluster evolution upon 
the various initial separation distributions used to define the 
VLMB population in the clusters. 

4-2.1 The logio -normal fit by Thies & Kroupa 

In one ensemble of clusters we select sep arations from the 
logio- normal fit to the VLMBA data by IThies fc Kroupal 
(2007) . In Fig. we show the results of dynamical evolu- 
tion in a dense cluster with a half-ma ss radius of 0.1 pc. 
When using the lThies fc Krou pa (2007) logio-normal sepa- 
ration distribution as an initial condition, i t is impossible to 
dynam ically disrupt many of the systems. IThies fc Kroupal 
(2007) used a binary fraction of /vlmb = 0.15 to fit the 
observational data. We show in Fig. [3] that the observa- 
tions can only be recovered if an initial binary fraction of 
/vlmb = 0.15 is used. Any other (higher) fraction hugely 
overpopulates the distribution around its peak, even when 
adopting a dense cluster model. 

This initial separation distribution does not account for 
the wide systems, which are observed in both young clus- 
ters and_jri_th«jieki. Nojnitial binary fraction coupled with 
the IThies fc Kroupal l|2007l ) logio-normal fit can reproduce 




0.01 0.1 



1 10 

Separation (au) 



Figure 3. The evolution of the separation distribution for 
VLMBs with initial se parations drawn from th e logio-normal fit 
to the observed data bv lThies fc Kro upaj l|20 07 | . the (solid) br own 
line). The wider fit to the data bv lBasriX: Reinersl l|2006l . the 
(dot-dashed) magenta line) is shown. The initial VLMB fraction 
is /vlmb = 0.15, and the half-mass radius of the cluster is 0.1 pc. 
The open histogram shows the initial distribution, and the hashed 
histogram shows the distribution after 1 Myr. For comparison , 
the logio-normal fits for field M-dwarfs l lFischer fc Marcvlll992i 
the ( dashed) blue line) and field G-dwarfs l|Duau^nn*ov**fcAlavo*r1 
Il99ll . the (dotted) red line) are also shown. 



them in any cluster . This is analogous to the result found by 
iKroupa fc Burkertl <|200lh . who demonstrated that it is im- 
possible to widen the stellar binary separation distribution 
through dynamical evolution in clusters. 



4-2.2 The Very Low Mass Binary Archive data 

In one subset of cluster models, we created the VLMB pop- 
ulation by selecting systems at random from the online 
VLMBA. This leads to the separation distribution shown 
in Fig. [D 

We dynamically evolve three ensembles of clusters con- 
taining these VLMBs. One is a very dense cluster with a 
half-mass radius of 0.1 pc, and an initial binary fraction 
/vlmb = 0.25. The second is a low-density cluster with 
a half-mass radius of 0.8 pc (similar to that of Orion to- 
day), and an initial /vlmb = 0.25; and the third is a low- 
density cluster with a half-mass radius of 0.8 pc, and an ini- 
tial /vlmb = 0.15. The results are shown in Figs. [4] and [5] 



In the dense cluster (ri/ 2 = 0.1 pc, Fig. 4(a) I, the wide 
VLMBs in the range 100 - 1000 au are all destroyed, and 
the very wide binaries (> 1000 au) cannot form at all. This 
is because the average distance between stars in clusters of 
this density is ~ 2000 au, making it highly unlikely that the 
wide binaries will be bound syste ms, even before dynamical 
evolution (see lParker et al1l2009l . for a detailed discussion). 
There is some disruption of the intermediate (4 - 100 au) 
VLMBs, but not enough to drastically alter the initial dis- 
tribution. 
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Figure 4. The evolution of the separation distribution for VLMB with initial separations and masses drawn from the VLMBA. The 
initial VLMB fraction is /vlmb = 0.25, and the half-mass radius of the cluster is (a) 0.1 pc and (b) 0.8pc. The open histogram shows 
the initial distribution , and the hashed histogram shows the distribution afte r 1 Myr. The log i p-nor mal fit to the observed data by 
iThies fc Kroupal l|2007l . the (solid) brown line) and the wider fit to the data bvlBas ri fc R einersl l|2006l . the (dot-dashed) magenta line) 
is shown. For comparison, the logio-normal fits for field M-dwarfs l lFischer fc Marc^ll99S the (dashed) blue line) and field G-dwarfs 
l|Duquennov fc M ayor 199l], the (dotted) red line) are also shown. 



When we adopt an initial half-mass radius of 0.8 pc 



(Fig. 4(b) I, all the VLMBs placed in the clusters are found 



to be bound systems by our algorithm. However, dynamical 
evolution acts to break up the widest (> 5000 au) VLMBs 
and systems with separations less than this are unaffected. 

It is clear that the simulations have too high an initial 
multiplicity if we are to match the observational data. In 
Fig. [5]we show the r 1 / 2 = 0.8 pc simulation with a lower 
initial multiplicity (/vlmb = 0.15). This produces the cor- 
rect distribution for the intermediate binaries, but under- 
produces the number of very wide (> 1000 au) binaries re- 
quired to be consistent with the observational data. 



4-2.3 The logio -normal fit by Basri & Reiners 

We also conduct simulations in which separations are cho- 
sen from the wider log ip-normal fit to the observed data by 
iBasri fc Reinersl l|200d ). in a dense cluster with a half-mass 
radius of 0.1 pc. The initial multiplicity is /vlmb = 0.25; we 
find that any higher multiplicity (e.g. /vlmb = 0.5) vastly 
overproduces the number of binaries compared to the obser- 
vations, even after dynamical evolution. 

The resultant separation distribution after 1 Myr for 
the cluster with /vlmb = 0.25 is shown in Fig. [6] There 
is some dynamical processing of binaries with separations 
in excess of around 1 au, so that in a dense cluster the 
I Basri fc Reinersl d2006l) separation di s tribut ion is not recov- 
ered. And, like the IThies fc Kroupal (|2007T ) separation dis- 
tribution, this distribution cannot produce the very wide 
systems initially (and even if it did they would be destroyed 
in a dense cluster). 
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Figure 5. The evolution of the separation distribution for 
VLMBs with initial separations and masses drawn from the 
VLMBA. The initial VLMB fraction is /vlmb = 0.15, and the 
half-mass radius of the cluster is 0.8 pc. The open histogram shows 
the initial distribution, and the hashed histogram shows the dis- 
tri bution after 1 Myr. The logio-normal fit to the observed data 
by IThi es & Kroupa (2007, the (solid) bro wn line) and the wider 
fit to the data bv lBasrifc Reinersl l |2006l . the (dot-dashed) ma- 
genta line) is s hown. For comparison, the logio-normal fits for 
field M-dwarfs jFisdier fc Marcvl Il992l, the (d ashed - ) blue line) 
and field G-dwarfs iDuauennovfc Mavorlll99ll . the (dotted) red 
line) are also shown. 
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Separation (au) Separation (au) 



(a) ri/2 = O.lpc (b) r 1/2 = 0.8 pc 

Figure 7. The evolution of the separation distribution for VLMBs with initial separations drawn from the same logio-normal distribution 
as the stellar binaries. The initial multiplicity of the VLMBs is /vlmBb = 0-5 and w e show the results for two d ifferent half-mass radii ; 
0.1 pc (a) and 0.8 pc (b). The logio-normal fits to the data bv lThies &: Kroupal J20071 . the (solid) br own line) and |Basri fc Reiner 
the (dot-dashed) magenta lin e) are shown. For compariso n, the logio-normal fits for field M-dwarfs l|Fischer fc Marcvll 19921 . the (dashed - ) 
blue line) and field G-dwarfs jPuquennov fc Mavor|[l99ll , the (dotted) red line) are also shown. 
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4-2.4 The logio-normal fit for field G-dwarfs 

In the final sub-set of cluster ensembles, the semi-major 
axes of the VLMBs are cho sen from the same distribution 
l|Duauennov fc Mavorl Il99lf ) as those of the stellar bina- 



ries. In Fig. 7(a) we show the separation distribution of the 



VLMBs initially (open histogram) , and after 1 Myr (hashed 
histogram) for a dense cluster (r±/ 2 = 0.1 pc) with an initial 
multiplicity of /vlmb = 0.5. 

The resultant dynamical processing of the binaries in 
this dense enviro nment causes th e sepa ration distribution 
to evolve into a iBasri fc Reinersl (|2006l ) logio-normal, al- 
beit with some over-production of binaries with separations 
< 1 au. 



In Fig. 7(b) we repeat the simulation, but adopt an 
initial half-mass radius of T\i% = 0.8 pc for the cluster. With 
this initial condition, the cluster is not as dense, and we 
overproduce VLMBs at all separations > 10 au as dynamical 
destruction is very ineffective in such low-density clusters. 



Figure 6. The evolution of the separation distribution for 
VLMBs with initial separati ons drawn from t he lo gip-normal 
fit to the observed data by [Basri & Rcincrs (2006, the (dot- 
dashed) magenta line). The initial multiplicity is /vlmb = 0.25 
and the initial half mass radius is O.lpc. The open histogram 
shows the initial distribution, and the hashed histogram shows 
the distributio n after 1 Myr. The log ip-normal fit to the ob- 
served data by lThie~ fc Kroupa (200?], the (solid) brown line) is 
shown. For comparison , the logio-normal fits for field M-dwarfs 
I Fischer_&_^larcy 1992, the (dashed) blue line) and field G-dwarfs 
I Duouennov^rMwOTlll99ll . the (dotted) red line) are also shown. 



4.3 Multiplicity fraction 

In each simulation, we determine the multiplicity fraction 
at each timestep. For simplicity (and owing to the very few 
higher-order systems that form in the simulations), we ig- 
nore triples, quadruples, etc. and simply determine the bi- 
nary fraction of the VLMBs. We show the evolution of the 
VLMB fraction for two of our simulations in Figs. [S]and[!|] 
When the initial binary fraction is 0.5 and s eparations 
are drawn from the iDuquennov fc Mavorl (|l99ll ) distribu- 
tion (Fig. |8]) , there is considerable break-up of (wide) sys- 
tems within the first few crossing times within a cluster with 
an initial half-mass radius of O.lpc (as detailed for stellar 
binaries m IParker et a"ill2009h . The initial binary fraction is 
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Figure 8. The evolution of the VLMB fraction over 10 Myr 
for a cluster with t he V LMB separations drawn from the 
iDuquennov fc Mavorl jl99ll) distribution and an initial binary 
fraction of /vlmb = 0.5. The cross shows the initial binary frac- 
tion. The initial half-mass radius of the cluster is 0.1 pc. 



Figure 9. The evolution of the VLMB fraction over 10 Myr for a 
cluster with the VLMB separations drawn from the VLMBA and 
an initial binary fraction of /vlmb = 0.25. The cross shows the 
initial binary fraction. The initial half-mass radius of the cluster 
is 0.1 pc. 



calculated to be 0.42, due to the widest binaries inputted 
into the simulations no t being bound within the dense clus- 
ter l|Parker et all 120091 ). The final binary fraction is 0.29, 
within the uncer t ainty associated with the determination of 
iBasri fc Reinersl |200rJ ) of 0.26 ± 0.10. If we assume a less 
dense initial cluster configuration of 0.8 pc (to enable the 
formation and preservation of the widest VLMBs observed 
in the field), then the final multiplicity becomes 0.4, which 
does not agree with the upper limit of the observed value. 

We also examine the change in binary fraction in a sim- 
ulation using the VLMBA data as initial conditions. During 
the cluster evolution, the initially lower binary fraction of 
/vlmb = 0.25 is also reduced (Fig. [9]), but only by 0.05. This 
is due to the majority of systems being close, and hence not 
susceptible to break up, but also a lower initial binary frac- 
tion will be dominated by the single very low mass objects 
- any break up of the small number of binaries will do little 
to the global binary fraction. 



4.4 The evolution of the VLMB mass ratio 
distribution 

In Fig. [TD]we show the effect of cluster evolution on the mass 
ratio distribution of the VLMBs. The results in the plot 
are for a cluster with an initial half mass radius of 0.1 pc, 
with separations and masses chosen from the VLMBA. The 
initial distribution is shown by the open histogram and the 
final distribution is shown by the hashed histogram. The 
distributions are both normalised to the total number of 
initial binaries. 

Fig. \TU\ clearly shows that even in a dense cluster, the 
dynamical interactions do not change the mass ratio distri- 
bution of VLMBs. 
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Figure 10. The evolution of the mass ratio distribution of 
VLMBs with initial masses and separations drawn from the 
VLMBA. The initial half- mass radius of the cluster is 0.1 pc and 
the data are normalised to the number of initial binaries. The 
open histogram shows the initial distribution, and the hashed 
histogram shows the final distribution. 



5 DISCUSSION 

We have examined the dynamical evolution of the binary 
fractions and separations of a separate VLMB population 
in star clusters. Our aim is to investigate what the initial 
VLMB multiplicity and separation distribution could have 
have been to reproduce the current observations. 

Before starting, there are two key points to be consid- 
ered in this discussion: 
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Firstly, many binaries are susc eptible to de s truction 
in clu st ers (as seen above, also see Heggid 19751; Kroupal 
Il995al lbl; iKroupa et all 1 19991 . 120031 : 1 Parker et alJl2009l ). How 
likely a binary is to be disrupted depends on its separation 
and the density of the cluster in which it was born. In low- 
density environments almost all VLMBs can survive, whilst 
in dense clusters VLMBs with separations > 100 au are sus- 
ceptible to disruption (see Fig. [7]). 

Secondly, the field is the sum of star formation from 
clusters of all masses and densities, and the field population 
has already been dynamically processed. Therefore, in the 
field we are not seeing the birth population of binaries, but a 
mixed and evolved population (|Goodwinll2010h . As the bulk 
of observed VLMBs are in the field we are not seeing their 
birth population. 



5.1 Separation distribution 

Several authors have attempted to parameterise the present- 
day separation distribution of the VLMBs. We set up the 
VLMBs in our clusters with four different initial separation 
distributions to investigate whether any of these evolve into 
the observed distribution, depending on the initial density of 
the cluster (defined by the half-mass radius) and multiplicity 
of the population. 

In the simulations in which we use the lThies &: Kroupal 
< |2007h logio-normal fit to the data as the initial VLMB sep- 
aration distribution, we find that dynamical processing in 
even the most dense clusters does not have a significant ef- 
fect on the initial binary population. The initial separation 
distribution does not, by construction, contain wide VLMBs, 
and these are not formed within clusters. However, this dis- 
tribution could be the initial VLMB distribution if the very 
wide binaries are able to form vi a capture during the cluster 
dissolution phase as proposed bv lKouwenhoven et al.l (|2010T ) 
for very wide stellar binaries. In this case the initial bina- 
rity of the VLMB population must be low at formation as 
few binaries are destroyed (and conversely a few very wide 
systems may be created). 

Using the currently observed VLMBA data as the ini- 
tial distribution also does not reproduce the observations. 
Again, the close binaries in the main peak are virtually un- 
affected by dynamical disruptions. In this case very wide 
VLMBs can survive in low-density clusters, but they are de- 
stroyed in high-density clusters. Therefore starting with the 
observed population cannot reproduce the observed popu- 
lation unless, once again, some very wide binaries are pro- 
duced during cluster dissolution. 

Adopting the much wider iBasri fc Reinersl l|2006h dis- 
tribution as an initial condition has some success in re- 
producing the observed population in that it starts with a 
number of wider VLMBs and so, even after dynamical pro- 
cessing some wider binaries remain. However, this distribu- 
tion requires us to believe that a large cl ose VLMB popula- 
tion exists which is currently unobserved dMaxted fc Jeffries! 



Joergens! l|200S(l 



l2005l ; iBurgasser et~ail l200?t ). though see 
and also that a significant population of VLMBs with sep- 
arations of tens of au also remains unobserved, which is 
rather more controversial as direct imaging surveys are sen- 
sitive to this separation regime in both clusters and the field 
l|Burgasser et alJl2007l : |joergensll200Sh . 

Using the current field initial M-dwarf separation dis- 



tribution as the initial VLMB distribution in part fails 
to match the current VLMB population. Whilst high- 
den sity clusters can be ev olved to something resembling 
the IBasri fc Reinersl (|200rJ ) distribution, low-density clus- 
ters hardly process their VLMB population at all resulting 
in far too many wide- and intermediate-separation VLMBs. 

Therefore we can conclude fairly robustly that the cur- 
rently observed VLMB fraction and separation distribution 
are not those of the birth population. This raises the obvi- 
ous question of what are the binary fraction and separation 
distribution of the birth VLMB population? To answer this 
we have to address two inter-related questions. 

Firstly, do all clusters produce the same birth popula- 
tions? And secondly, what is the mixture of cluster densities 
that contribute to the field? It is the combination of cluster 
densities rather than masses that is important, as it is den- 
sity that controls dynamical processing. Any high-density 
cluster will process its wide binaries, whilst a low-densi ty 
cluster will not, regardless of its mass (|Parker et al.ll2009h . 

If all clusters do not produce broadly similar birth pop- 
ulations then it becomes almost impossible to draw many 
conclusions about birth populations other than that after 
dynamical processing they sum to make the currently ob- 
served population (see Goodwin 2010). However, if we ac- 
cept as a starting point that star formation is roughly uni- 
versal, and that the birth populations of VLMBs in all clus- 
ters are roughly the same then we are able to draw some 
conclusions. 

The pertinent points are all illustrated in Fig. [7] In 
this figure we see two key points: high-density clusters are 
extremely efficient at destroying VLMBs with separations 
> 20 — 100 au, whilst low-density clusters leave their birth 
populations relatively intact; and that no cluster can effec- 
tively alter the birth populations of VLMBs with separations 

< 20 au. Therefore, the VLMB population with separations 

< 20 au must be very close to the birth population, and the 
birth population of VLMBs with separations > 20 au will 
have been dynamically processed and must have contained 
many more systems initially than are observed in the field. 

This implies that the (putative universal) birth VLMB 
fraction is higher than we now observe, and that the birth 
separation distribution must ha ve a significant excess of sys- 
tems with wide separations (see lKroupalll995bl . for a similar 
argument for stellar binaries). 

How much higher the birth binary fraction must have 
been depends on the density distribution of star clusters. 
It should be noted that this is not the current density dis- 
tribution, but rather the range of maximum densities clus- 
ters have reached during their evolution. There are sugges- 
tions that many clusters undergo an early dense phase in 
which their densities can reach those simulate d by the 0.1 pc 



half-mass radius clusters we model here (se e Kroup a et al 
1999;lMoraux et all 2007; Bastian et ai]|200Sl; lAllison et al 



2009; iParker et alTl2009l ; iGoodwinl 120101 ). If many clusters 
do indeed undergo an early dense phase, then the initial 
VLMB population may need to be similar to the very wide 
/mult = 0.5 population illustrated in Fig. 7(a) If, on the 



other hand, many stars are created in low-density associa- 
tions, the initial VLMB population would only need to con- 
tain slightly more wide binaries than are now observed. 

What is clear is tha t VLMBs and M-dwarfs (taking the 
iFischer fc Marcvl (|l992r i distribution) must have had very 
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different birth populations. M-dwarfs are subject to exactly 
the same dynamical processing as VLMBs (only slightly 
modified by their larger masses), and to have such signif- 
icantly different separation distributions in the field means 
they must have been born with significantly different distri- 
butions. 

The key regime in which more observations are required 
is in examining if there are more wide and, especially, in- 
termediate separation VLMBs. The confirmation of either 
presence or lack of a significant VLMB population with sep- 
arations of ~ 100 au where we know that there are significant 
numbers of M-dwarf binaries is crucial in order to determine 
how different the initial binary populations are. 



5.2 Mass ratio distribution 

As demonstrated in Fig. 1101 the effect of dynamical evolution 
on the mass ratio distribution of the VLMBs is negligible, 
even in the most dense clusters. What is clear is that there 
is no preferential break-up of systems with a particular q 
value; the distribution is uniformly lowered by dynamical 
processing, in agreement with the evolution of the mass ratio 
distribution for low-mass stars (0.1 Mq 5j m p =^1.1 M^ and 
0.1 Mq ^ m s ^ 1.1 M©) in clusters found by iKroupa et al.l 
(|2003h . This means that the current mass ratio distribution 
represents the birth mass ratio distribution. 



6 CONCLUSIONS 

We present the results of iV-body simulations of the effect of 
dynamical evolution on very low mass binaries (VLMBs) in 
'typical' N — 2000 Orion-like star clusters with various ini- 
tial densities. In each cluster, we place a separate population 
of VLMBs and follow the effects of dynamical interactions 
on these systems. 

Our conclusions can be summarised as follows: 

• VLMBs with separations < 10 — 20 au cannot be dis- 
rupted, even in the densest clusters. Therefore the currently 
observed VLMB population with separations < 10 — 20 au 
must be the primordial population. However, we note that 
the VLMB population at these separations is poorly known. 

• Many VLMBs with separations > 100 au can be 
destroyed in very dense clusters (with half-mass radii of 
0.1 pc), but are relatively unaffected in low-density clusters 
(with half-mass radii of 0.8 pc, so 512 times less dense). 

• If all star clusters produce the same birth VLMB pop- 
ulation, then the birth VLMB binary fraction must have 
been somewhat higher than we observe in the field. The ini- 
tial VLMB distribution must also have significantly more 
wide VLMBs than are currently observed unless capture 
during cluster dissolution is an effective wide binary for- 
mation mechanism. 

• The mass ratio distribution of VLMBs is unchanged 
by dynamical processing and so is a probe of the birth mass 
ratio distribution. 

• M-dwarf binaries are also processed by clusters in the 
same way as VLMBs, and so the very different field popu- 
lations must reflect very different birth populations. 

Further work to better constrain (and improve the 
statistics of) the separation distribution of both VLMBs and 



M-dwarf binaries in the field would be beneficial in determin- 
ing whether the transition in separations between M-dwarf 
and VLMBs is as abrupt as Fig. [1] suggests. If not, then 
the trend in decreasing multiplicity with primary mass, and 
fewer wide VLMBs than wide G-dwarf binaries, suggests a 
possible continuation through the hydrogen-burning limit 
from the stellar to the substellar regime. 
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